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ANNALS OF M ATHEMATICS. 

Vol. III. December, 1887. No. 6. 

ON THE ORBIT OF HYPERION. 
By Prof. Ormond Stone, University of Virginia, Va. 

1. The principal difficulty in the integration of the equations of motion in 
the case of the problem of three bodies arises in the integration of terms involv- 
ing the inverse powers of the distance between the disturbed and disturbing bod- 
ies. When the ratio between the radius vectors is not too near unity, the inverse 
powers referred to can be developed in rapidly converging series, in terms of 
multiples of the elongation. When, however, these ratios do not differ greatly 
from unity, the convergence of the series mentioned is very slow. If, in addition, 
the mean motions of the two bodies are nearly commensurate, the ordinary 
methods of solving the problem become inapplicable. 

Such a case presents itself in the determination of the perturbations of Hype- 
rion produced by Titan. On the other hand, the mutual inclination of the orbits 
of these two satellites is so small as to have eluded detection ; the eccentricity of 
the orbit of Titan is less than 0.03 ; and the position of the apo-saturnium of Hy- 
perion so nearly coincides (at least at present) with the point of conjunction of the 
two satellites, as to give rise to a suspicion that the eccentricity of its orbit is 
small, and that the apparent eccentricity is principally due to the perturbations 
produced by Titan. 

2. First Approximation. — These circumstances have suggested* the propriety 
of first investigating that part of the disturbance which may be determined by 
neglecting the mutual inclination and the eccentricities of the orbits of both bodies. 
On this hypothesis Tisserand has found that, assuming the mass of Titan 
m' = ^ft^gfl , the motion of Hyperion may be represented, to terms of the first 
order with respect to the mass, by the formulai, 

- = I — 0.0004 cos (/' — /) — 0.0014 cos 2 (/' — /) 
+ 0. 1000 cos 3 (/' — /) + 0.0006 cos 4 (/' — /) , 

♦Tisserand. Sur un Cas Remarquable du ProbUme des Perturbations. Comptes-Rendus, 1886, 2» 
Semestre, p. 446. 

Hill. Coplanar Motion of Two Planets, One Having a Zero Mass. Annals of Mathematics. 
Vol. in. p. 65. 
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w = /-{- lo'sin (/'—/)+ 13' sin 2 (/' —/) 
-683' sin 3 (/'-/)- 3' sin 4 (/'-/), 
where r is the radius vector of its orbit, 

a mean distance from Saturn, 

w true longitude in orbit, 

/ mean longitude in orbit, 

/' mean longitude in orbit of Titan. 

3. Second Approximation. — As a second approximation let us assume 

r^a(\ -\- a^ cos 9 -\- a^ cos 26 -\- . . .')^= a(\ -\- a), (2) 

— = «(i + «, cos ^ + ;«2Cos 2 ^ + . . .)^«(i+r), (3) 

where r and w are the radius vector and longitude in orbit of Hyperion, 
d ^P — l=(n' — «) / is the mean angular distance between the radius vectors 
of Titan, and Hyperion, a and « are the mean values of r and w, and a,, «i, etc. are 
constants to be determined. 

The differential equations of motion may be written 

,dw 



d}- 
dt 



r* ^ = m'J^jSrdt -f C, (4) 



where k is the constant of the system as derived from the orbit of Titan, m' the 
mass of Titan, tn'k^R and m'l^S the components in the plane of the orbit of the 
disturbing force in the direction of and perpendicular to the disturbed radius vec- 
tor, and C is a constant of integration. As is well known, when the orbit is 
nearly circular, Cis approximately equal to k\/a\ say k\/\a{\ — i*)], where v is 
a small quantity. 

4. Substituting in (4) the values of r and dw\dt given by (2) and (3), and di- 
viding by (^n, we have 

r^ dw , , , . , ,2 k\/{\ — v) , , k^ C <r J^ /a.\ 

■^^-^=i+2o + r + d^-^2ar + ah= ^ ^,^ ) + ,n' ^J Srdt, (6) 

in which, if 2" = I, i being integral and positive, 

<r = rtj cos d ^ a^ cos 2d -\- . . . = Zai cos iD, 
7 = n^ cos d -\- n, cos 2^ + . . . = 2'«j cos it), 
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<t' = ^2'ai^ -t ^'«irti + 1 cos d + (^«i' + I'atai+i) cos 2d 

+ (a,n^ + i'«i«i + 3) cos 3# + (^«2* + a.a, + i'«i«i + ,) cos 46 
+ (rt,«4 + ajiJj + I'^i^j^..,) COS 5<^ + . ... 

5. Comparing (2) and (3) with Tisserand's values of rf a and w, we are led 
to infer that it will be proper, for the present, to consider a, and «, as quantities 
of the first order ; a^. a^, «<, «„ »,, n^, as quantities of the third order; «„ a^, cl,, a^, 
«5, «„ «j, fig as quantities of the fourth order ; and neglect the remaining coeffi- 
cients of a and r. 

We have, then, to terms of the fourth order, 

<T ^ <?i cos ^ + «2 cos 26 -\- . . . + tfg cos 8^, (7) 

r = «, cos ^ + «2 cos 2I) -\- . . . + «g cos 8<?, 
<j2 =: ^a^ 4- aij (d-j + a^ cos ^ + a,rt, cos 2d + aia^ cos 46 (8) 

4- «2«3 cos $0 + ^«/ cos 6^ + a^a^ cos 7^, 

2<r7- = rt3«j + [(«2 + '^4) «3 + «3 («j + «<)] cos (? 

+ (^J^i^s + ^t^h) cos 2^ + (ai«3 + rtjW,) cos 4^ 
+ (^j'^a + <^3«2) cos S^' + aj«3 cos 6^ 

+ {«3«4 + «4«3) COS 7^, 

a'r = ^a^n^ cos 3^ + \a^n^ cos 9^. 
Whence (6) becomes 
r* dw ,12, 

+ [2«i + «i + rta («2 + «4 + «2 + «j) + (^2 + «4) «3] COS 
+ [2^2 + «2 + «S («l + «l) + "^I^s] COS 2d 

-\- [2^, + »3 + \a^n^ cos 3(9 

+ [2^4 + «4 + «3 («, + «i) + rt,«3] cos 4<^ 
+ [2rt5 + «5 -J- «3 (^2 + «2) + ''2«3] cos 5</ 
-j- [2«j + «6 + \(^i + '^3«3] COS 6d 

+ [2^7 + «7 + «3 (^4 + «<) + a^n^ cos 7(y 
-j- [2rtg + ^g] cos 8^ 



^fn'^JSrdt. (9) 
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6. Remembering that C= k}/[^a (i — v)], equation (4) also gives 

whence (5) becomes -rfl + ^\r — a{\ — v)\ = m'/^P, 

where P = R -\- p\^ Vi" {1 — v)] . 2^Srdt + m'/^ f f^n^/l '] ; (10) 

orsmcer = a{i+a), ^ + ^^^—^ = „,' - R („) 

7. Differentiating twice with regard to t, since 6 = {n' — n) t, where «' is 
the mean motion of Titan, this gives 

-y-j- = — («' — nf («, cos ^ + 4 «2 cos 26 + 9«3 cos 3^ + . . . ) . (12) 

Expanding by the binomial theorem, 

*^ + " ; = V + (i — 3v) <r — 3 (i — 2v) tr" + 2 (3 — 5v) ^r' — 5 (2 — 3v) ff* 



+ 3 (5 - 7»') '^ - 7 (3 - 4^-) '^ + • • • • 

In the development of the third and higher powers of <r, it will be sufficient 
to put <T = «3 cos 3^ ; whence 

(j^ = \a^ cos 1(1 + \ai cos 9^, 

<7* = ^ ^3* -)- \a^ cos 6^ + ^^3* cos 1 2d, 

(13) 
(T* = \a^ cos 3(? + 35^a3= cos 9^ + ^^3' cos 15^, 

<T^ = ^«3^ + ^a^ cos 6^ + A^s' <^os 12^ + ^^"''^a' cos i8#. 

8. If we neglect terms whose arguments are multiples of ^ greater than 8^, 
we have to terms of the fourth order 

"iX^of = ^' — 1^3' (I - 2v) - ^-a^' (2 - 3^) 
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+ [«1 (l — 3*') — 3«3 («2 + '^O (l — 2V)] COS 

+ [«2 (i — 3v) — 3a,«j (i — 21-)] cos 2(? 

+ i'^i (i — 30 + W (3 — S"^)] cos 3d 

+ [«4 (i — 3*') — S'^i^s (i — 2v)] COS 4<? (14) 

+ K (i — 3*') — 3«2«3 (i — 2v)] COS id 

+ K (i — 31-) — W (i — 2v) — f^a* (2 — 3v)] COS 6^ 

+ [«, (i — 3v) — 3«3«< (i — 2v)] cos 7(? 

+ «8(i — 3v) cos 8^. 

Since the left hand members of (9) and (11) contain only terms involving 
■cosines of multiples of and constant coefficients it is evident that the right 
hand members of the same contain similar terms only, and we may put 



. , , ; i Srdd = 2 Si cos id, 

(tn («' — n)J , 



o 
•whence equations (9) and (11) may be written in the form 

9. With assumed values of our constants, R and .S" are computed for partic- 
ular values of 0, selected at equal intervals from ^ ^ o to ^ = 180°, and then ex- 
panded mechanically into a series of cosines of multiples of 0, effected with con- 
stant coefficients ; whence F and j Srd0, and thence Pi and Si are easily obtained. 

If now we replace the left hand members of equations (15) and (16) by the 
developed expressions of (9) and (11), and equate the coefficients of the cosines 
of equal multiples of 0, there will result any number of equations from which 
may be derived a corresponding number of the quantities «j, w,, etc. Instead, 
however, of considering the mass of Titan as known, it will be better to assume 
^3 as given, and to consider m' as one of the unknowns. 

With Tisserand's values of the coefficients, the constant part of (16) gives an 
approximate value of v; with this value of v the constant part of (15) gives an 
approximate value oi fi = c? {n' — nf / P ; with these values of [i and v, the terms 
of (16) involving 3^ give an approximate value of;;/'. On account of the impor- 



'^ y={"^^-]\i-^). (18) 

h»4 + i«3X j ^ » J 
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tance of these quantities I have developed the equations giving u, m' , and v to 
terms of the sixth, seventh, and eighth orders respectively. The remaining equa- 
tions have been carried only to terms of the fourth order. 

ID. The equations giving v, //, and m' are as follows, including terms of the 
orders named : — 

V [i + 3 («.' + <^2 + 'j/ + «4*) + ¥«3* + ^a,* + \'Wi^ (17) 

f I + i («,^ + «/ + «/ + a,^) 

( + «1«1 + «2«2 + «3«3 + ^i' 

m'P^ = (i — 3v — 9/^) ^3 _ 3 (i _ 2v)>, («j + a^) 

+ t (I - f v) «/ + I^ (I - |v) «/ + 2^ (I _ ^v) a/. ^'^^ 

Comparing (16) with (12) and (14) the coefficients of cos 6, cos zd, cos 4^ 
give the following equations for the determination of a^, a^, a^ : — 

(1 — 3v — /^) «, — 3*, (i — ai-) «j — 3^3 (i — 2v) a^ = m >P^, 

— 3^3(1 — 2v) «, + (i — 3v — 4/i) a, = w'Pj, (20) 

— 3^3 (i — 2v) a, + (i — 3" — i6/i) a^ = m '/"<. 
Similarly, for a^, . . . a^, we have 

(i — 3v — 25//) «5 = :;«'/'5 + 3<?s (i — 2V) «j, 

(i — 3v — 36//) «5 = Wi'Ps + 3^3 (i — 2v) . ^^3 + 5«3*(i — |v), 

(21) 
(i — 3v — 49/i) «, = w 'P, + S'^s (i — 2v) rt^, 

(i — 3v — 64/i)t?, = w'/'j. 

In the same way «3 is obtained from 

(i + i'^i) «3 — m'Si — 2at, (22) 

which is derived by comparing (15) with (9); after which «,, «,, and n^ are ob- 
tained from 

«1 + rt3«j + a3«i = w'^i — 2rt, — (ff, + Ws) {a J + ilj, 

«73«i + Wj =m'Si — 2^2 — (ffs + Ws) rti, (23) 

iZjW, -(- n^ = m' St — 2rt< — (rt, + n^) a, ; 
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and the remaining terms give 

«, = tn'S^ — 2^5 — (a, + «s) a^ — aiti^, 



n^ = m'S^ 



2a. 



W — «s«s, 



«7 = »t'S, — 2aj — rtj (a^ + «<) , 



(24) 



»g = m 'Sg — 2ag. 



With the values of the constants thus obtained, new values are determined. 
This process must be repeated until the values obtained agree with those assumed. 
1 1. For the determination of R and 5 were employed the following : — 

n = i6?9i99, n' = 22?5770, 

- = I — 0.0004 cos d — 0.0014 cos 2d 

+ o.iooo COS 3^ + 0.0006 cos 4<?, 
w' — w = d — 10' sin d — 13' sin 26 
+ 683' .sin 3(? + 3' sin 46, 
log a = 0.0835122, r' = 1, 

f? = >-'' sin^ iw' — w) + \r' cos iw' — w) — r^, 

t —L _L 
« — ^s ^/s . 

y 

R = hr' cos (w' — w) 3 , 

5 = hr* sin (zc' — w) . 

With these data values of a^R and Sr were computed for every 20° from 
^ = o to (? = 180°, as follows : 



« 




a^R 


5/- 


0° 


— 


14.815 


1 
0. 


20 


— 


S-.U/ 


+ 2.170 


40 


— 


1.926 


+ 0.340 


60 


—_ 


1.501 


— 0.134 


80 


— 


1. 1 30 


— 0.623 


100 




0.438 


— 1. 182 


120 


+ 


0.414 


— 1-231 


140 


+ 


0.97s 


— 0.699 i 


160 


+ 


1. 126 


— 0.220 


180 


+ 


II33 


0. i 

■ 



p 
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Whence 

a'R = — 1.423 — 3.480 cos d — 1.538 cos 2d — 1.466 cos 3/? 

— 1. 551 cos 4^ — 1.322 COS 5^ — 1. 213 cos 6^ 

— 1. 149 cos yd — I.I 16 cos 8^ — 0.553 cos qd^ 
Sr= — 0.561 sin d + 0.800 sin 20 + 0.654 sin $9 

+ 0.3 II sin 46 + 0.459 sin 5^ + 0.368 sin 6d 
+ 0.232 sin 7^ + 0.131 sin 8^, 
SrdO = + 0.561 cos d — 0.400 cos 2^ — 0.218 cos 3^ 

— 0.078 cos 46 — 0.092 cos 5^ — 0.061 cos 66 

— 0.033 cos yd — 0.016 cos 8^. 

Substituting in aV"^ = i — 3<t + 6*7^ — lOfT* + . . , 

the values of <r, a^, etc., derived from (7), (8), and (13), employing Tisserand's val- 
ues of ax, a^, «s, etc., we have 

a^r^^ =^ 1.031 + 0.00 1 cos 9 + 0.004 cos 26 — 0.308 cos 3^ 

— 0.002 cos 4^ — 0.00 1 cos 5^ + 0.03 1 cos 69 

— 0.003 cos 9^. 
Whence (10) gives 

a*P=^ — 1. 241 + 0.353 cos — 4.399 cos 29 — 2.75 1 cos 7,9 

— 2.545 cos 49 — 1.457 COS 5^ — 1.386 cos 69 

— 1.224 cos 7^ — 1. 172 cos 8^ — 0.5 18 cos 9^. 
Differentiating (i) and comparing with (3), we have 

«i ^ + o.ooio, 
«j = + 0.0025, 

«3= —0.1993, 

n^= —0.0012. 

From equations (17) to (24) were next obtained the values given in the first 
column of the following table. With these and the values of P< and ^^ already 
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employed new values were obtained, which are given in the second column of the 
table: 





I 


II 


V 


0.01493 I 


0.01409 1 


log // 


9.054927 


9-055352 


ijm' 


1290. 


1370 


«i 


— 0.0013 


— 0.0012 


«2 


0.0076 


— 0.007 1 


«3 


+ O.IOOO 


+ O.IOOO 


«* 


+ 0.0027 


-f 0.0025 


«5 


+ 0.0018 


+ 0.0017 


a» 


+ 0.0003 


+ 0.0003 


«r 


0.0 


0.0 


«8 


+ 0.0002 


+ 0.0002 


«1 


+ 0.0027 


+ 0.0025 


»t 


+ 0.0137 


+ 0.0129 


«3 


— 0.1990 


— 0.1990 


n^ 


— 0.0062 


— 0.0058 


«5 


— 0.0023 


— 0.0022 


«6 


— 0.0008 


— 0.0008 


»7 


+ 0.0002 


+ 0.0002 


»8 


— 0.0003 


— 0.0003 



The corresponding formulae for the radius vector and longitude of Hyperion 



are 



- ^ I — O.OOI2 COS — 0.0071 COS 2d + O.IOOO COS 3^ 

+ 0.0025 COS 4/) -f 0.0017 COS 5*9 + 0.0003 COS 6d (^5) 

+ 0.0002 COS Sff, 
Zf = / -f- 26' sin ^ + 66' sin 26 — 682' sin 3^ 
— 1 5' sin 4^— 5' sin 5^ — i'sin6^. 



(26) 



I/O STONE. ON THE ORBIT OF HYPERION. 

Finally from the formulae 

and the values of ,« and in' obtained from the second column of the above table, 
we have 

log k = 1.353507, log a = 0.085727. 

12. As a check, the computations were carried through a second time, this 
time dividing the circle into 72 parts, and working independently of Tisserand's 
values of the smaller coefficients in the expressions for rj a and dw\dt. In com- 
puting the forces, I have assumed 

r =^ a(\ + o. I cos 3<^), 

w' — w = <^ -f 683' sin Tfi. 

To find a, if we neglect small terms of the sixth order and those containing 
;«', we have 

(28) 



2«3 



(?~n'^- i-v ' ^^^' 

whence «3 = 0.1485, v ^ 0.01493, 

a' = 1. 00000, log a = 0.08569. 
Including only those terms which will be needed in computing m', we have 
a^R^ — 1.280 — 1. 141 cos 3^ — . . ., 
I SrM = — 0.206 cos 3^ — 0.083 cos 60 — 0.038 cos gd — . . . , 
a*r~' ^ 1.031 — 0.308 cos 3^ + 0.031 cos 66 — . . . j 
0^= — 1.061 — 2.600 cos 3^ — . . , . 
Finally, the equation 

[3«3/'o + (I + ias)^^] ^' ^ [^7(7-!r^) -9/^(1+ W)J«3 , (3 1) 

obtained by substituting (17) in (19) and neglecting small terms of the sixth 
order, gave 
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13. This approximate value for the mass of Titan is about eight times 
that given by Tisserand. The reasons for this are apparent : To terms of the 
first order, n = («' — nf\n^ and m' P^^^ii — 9//) a^, P^ being obtained from the 
hypothesis r = a and r' = a'. tn'P^, however, is so small that the introduction 
of a more accurate value of ji and the addition of the remaining terms given in 
equation (19) change it {in'P^) from — 0.0006 to — 0.0021, while P^ changes from 

— 6.54 to — 2.60. The disturbance produced by Titan is, of course, greatest 
near conjunction ; but Hyperion is then always at apo-saturnium, and accord- 
ingly the average distance between the two satellites is approximately a{l -f ^3) 

— a' instead oi a ■ — a' ; i. e. to terms of the first order, 0.275a instead of o. 17512. 

14. Prof Newcomb in a discussion of the motion of the peri-saturnium of 
Hyperion*, assuming, as I have done, that the orbit of Titan is a circle, has 
found m' = xy^^cg- ; but in so doing has assumed n to be the elliptic mean mo- 
tion, instead of the mean motion in longitude: whereas the mean value of the 
mean motion in the variable ellipse is, approximately, kl afi. Assuming 

k 
k = n' and n ^ -j ^ i6°.793i, 

we have 4n — 3«' = — o°.5585 ; 

whence, using the same values for the forces as those employed by Newcomb, 

"' ~ 40.6 n -T^' 
which agrees approximately with the values already found. 

*On the Motion of Hyperion. A New Case in Celestial Mechanics, pp. 365 7. 



£tO BE CONTINUED.] 



